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ABSTRACT 

The two reaction wheel K2 mission promises and has delivered new discoveries in the stellar and 
exoplanet fields. However, due to the loss of accurate pointing, it also brings new challenges for the 
data reduction processes. In this paper, we describe a new reduction pipeline for extracting high 
precision photometry from the K2 dataset, and present public light curves for the K2 Campaign 1 
target pixel dataset. Key to our reduction is the derivation of global astrometric solutions from the 
target stamps, from which accurate centroids are passed on for high precision photometry extraction. 
We extract target light curves for sources from a combined UCAC4 and EPIC catalogue - this includes 
not only primary targets of the K2 campaign 1, but also any other stars that happen to fall on the 
pixel stamps. We provide the raw light curves, and the products of various detrending processes aimed 
at removing different types of systematics. Our astrometric solutions achieve a median residual of 
~ 0.127". For bright stars, our best 6.5 hour precision for raw light curves is ~ 20 parts per million 
(ppm). For our detrended light curves, the best 6.5 hour precision achieved is ~ 15 ppm. We show 
that our detrended light curves have fewer systematic effects (or trends, or red-noise) than light curves 
produced by other groups from the same observations. Example light curves of transiting planets and 
a Cepheid variable candidate, are also presented. We make all light curves public, including the raw 
and de-trended photometry, at http://k2.hatsurveys.org. 

Subject headings: K2, astrometry, photometry 

1. INTRODUCTION 

The Kepler spacecraft ended its primary mission after 
the failure of two reaction wheels. The K2 mission uses 
the Kepler spacecraft to perform 80-day observations of 
selected fields in the ecliptic plane. This brings new op¬ 
portunities to study transiting planets around different 
stellar populations compared to the original Kepler field, 
such as clusters of young and pre-main sequence stars 
(jHowell et al.ll2014l) . 

K2 uses the remaining two reaction wheels, and solar 
radiation pressure, to maintain close to constant pointing 
of the spacecraft over the 80-day per-field observations. 

Currently, observations are performed with 21 modules, 
each module consisting of 4 CCD channels, yielding 76 
channels (2 modules failed). Due to the limited band¬ 
width, only postage stamps containing proposed targets 
are downloaded. These postage stamps are typically 
25 x 25 pixels in size (depending on the brightness of the 
targets and campaigns). These make up only less than 
10% percent of the entire field of view (FOV). The ma¬ 
jority of stamps are observed at ~ 30 minutes cadence. 

Typically, two Full Field Images (FFIs) are downloaded 
for the beginning and the end of campaign. 

However, the two reaction wheel mode also brings in 
new challenges for the data reduction processes. The 
spacecraft pointing is less stable compared to the primary 
mission, leading to a potential decrease in the photomet¬ 
ric precision. Although the disturbance from the solar 
pressure is mostly controlled by the two reaction wheels 
and the thruster firing (every 2 days), there is still a 
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low frequency motion remaining, resulting in the targets 
drifting across the field of view. The extracted aperture 
photometry light curves are dom inated by t h e syst e mat- 
ics in duced by this drift pattern. iVanderburg fo Johnsonl 
(f20ll (hereafter VA14) minimized this drift system¬ 
atic by decorrelating the light curves with the motion 
of the spacecraft. They achieved a photometric preci¬ 
sion that is within a factor of two of the original Kepler 
photometry. Various other teams also develope d their 
own tools to reduce the K2 data. lAigrain et al. ( 12015t l 
used aperture photometry and a semi-parametric Gaus¬ 
sian process mode l to extract p hotom etry from the K2 
engineering data. iLund et all (I2015T ) presented K2P, a 
pipeline specifically de s igned for a strometric analyses . 
iForeman-Mackev et all 1 2015 1 and lAngus et all 6201 5(1 
proposed a method to analysis the K2 data without a 
general detrending process. 

There is, however, room for further improvements. 
VA14 reduction achieved the highest precision among all 
the past works, but only derived photometry for the pro¬ 
posed Kepler targets (not all targets falling on silicon), 
and are also known to have remainin g system atic varia ¬ 
tions affected by the spa cecraft roll (Angus et akl f20151. 
lAigrain et al.1 (|2015fl and ILund et al (1201511 derive pho¬ 
tometry for all of the targets on silicon, but achieved 
slightly lower precision than VA14, especially for the 
bright stars. Here we present a new reduction of the 
K2 data drawing on techniq ues used in analysing data 
from ground-based surveys dBakos et al.ll20ld e.g.). 

We approach the K2 pixel file reduction with the fol¬ 
lowing steps: 1) improved astrometry for source centroid- 
ing and flux extraction; 2) photometric extraction for all 
the stars observed on the K2 postage stamps; 3) removal 
of first order systematics via a modified External Pa¬ 
rameter Decorrelation (EPD) procedure (broadly similar 
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to VA14); 4) further reduction of the shared systematic 
trends via an implementation of the Trend Filtering Al¬ 
gorithm (TFA) and semi-periodic stellar oscillations via 
cosine-filtering. The global astrometry step is key to this 
process - it minimizes the effect of spacecraft drift on the 
aperture photometry, and allows us to accurately model 
the spacecraft motion for further detrending. 

In this paper, we describe our K2 photometry pipeline 
and the high precision light curves from the reduction of 
K2 Campaign 1. We introduce our effort of deriving ac¬ 
curate astrometry for the K2 observations, making use of 
the K2 FFIs, and present a revised K2 Campaign 1 tar¬ 
get list in §2. In §3, we present our aperture photometry 
method. In §4, we revisit our detrending techniques and 
present our light curves at different detrending stages. 
In §5, we compare our photometry with that of other 
studies. 

2. ASTROMETRY 

2.1. Background 

The first step of our reduction is to derive an accurate 
astrometric solution of the K2 data. Despite its large 
pixel scale (~ 4") and PSF FWHM (5 — 6 "), the origi¬ 
nal Kepler mission turned out to be a great tool for accu¬ 
rate astrometry itself because of its extremely high SNR 
photometry and stable pointing. iMonet et all (|2010h re¬ 
ported a preliminary astrometric solution precision, from 
the first few months of Kepler data, to be 0.001 pixel, 
nominal 4 mas. This high astrometric precision, and 
high stability of the centroid position, enabled the high 
photometric precision of the Kepler primary mission. 

Unlike the original Kepler Mission, the K2 stars typi¬ 
cally drift across the CCD plane at a speed of 1-3% of a 
pixel every 30 min. Since the 30 min K2 frame is com¬ 
posed of 270 short exposures of 6 s each, the final PSF 
is inevitably distorted, and neighbouring stars tend to 
become blended. Therefore, it is difficult to determine 
accurate centroids from source extraction alone. Thus, 
we use an external catalogue, namely the fourth United 
States Naval Obs ervatory (USNO) CCD Astrograph Cat¬ 
alogue, UCAC4 (jZacharias et, al.l[2013fl . which has an as¬ 
trometric precision of 15-100 mas, to derive good astro- 
metric solutions for the K2 frames. 

A good astrometric solution does not only benefit the 
photometric precision, but also enable us to make max¬ 
imal use of the K2 observations. The K2 campaigns ob¬ 
serve targets proposed by the community, and each tar¬ 
get was then assigned a stamp of size 20-50 pixels across. 
This stamp size is much bigger than the original Kepler 
stamp size. In addition to the target, many other sources 
are observed in a typical K2 stamp. We provide position 
information and reduced light curves for all of the stars 
observed in the K2 stamps. We anticipate an improved 
planet yield from this approach, due to the larger num¬ 
ber of sources available, and the availability of the light 
curves of neighbouring stars, useful in blend analyses. 

We first derive a general astrometric solution using the 
Kepler FFIs and our custom developed astrometry soft¬ 
ware used for HATNet. This general solution is then used 
as an initial guess for the remaining stamp observations. 
We stitch all the K2 stamps together into a “Sparse FFI” 
(SFFI), with the unobserved regions masked. We fit for 
an astrometric solution to the SFFI, which is assumed 


to be a low order polynomial distortion from the FFI 
astrometric solution. 

2.2. Astrometry Standard Catalogue 

We use the UCAC4 catalogue (|Zacharias et a. 1.1 [20 1 31) 
as our astrometry standard for deriving the astromet¬ 
ric solution. It contains over 113 million objects, and is 
complete down to magnitude R = 16. 

The precision of coordinates provided by UCAC4 is 
~ 15 — 100 mas. UCAC4 catalogue also contains proper 
motion of ~ 105 million stars, with errors around 1 to 10 
rnas/yr. Both the coordinates and proper motions are 
measured on the International Celestial Reference Sys¬ 
tem (ICRS) at a mean epoch of 2000. We linearly cor¬ 
rected the coordinates based on the proper motions to 
epoch 2014. UCAC4 also contains Two Micron All-Sky 
Survey (Y2MASS: iSkrutskie et al.lf2006l) ') photometry for 
around 110 million stars, and AAVSO Photometric All- 
Sky Survey (APASS) five-band (BVgri) photometry for 
over 51 million stars. For the stars in 2MASS but with¬ 
out APASS photometry, their gri band photometry are 
estimated using the 2MASS magnitudes. BV band mag¬ 
nitudes are adopted from the Tycho-2 catalogue where 
available, otherwise also estimated from 2MASS. We use 
the B and V magnitude to estimate the magnitude of 
stars in the Kepler band when needed. 

2.3. Astrometry on the Full Frame Image 

The Kepler FFIs are divided into subimages by read¬ 
out channels. There are 84 subimages for each FFI. Two 
of the CCD modules (8 channels altogether) failed dur¬ 
ing the Kepler main mission. The remaining 76 subim¬ 
ages were used to create the images from the 38 working 
CCDs (follo wing the Kepler Instrument Handbook). We 
use FISTAR (jPdl l2012ll for source extraction. The un¬ 
certainties of the source extractor is about 0.07 pixels. 
This is estimated by comparing the extracted source po¬ 
sitions on the two different FFIs taken from the Cam¬ 
paign 1. The relative shift and rotation between the 
two FFIs were taken into account by fitting a low order 
polynomial to the two extracted source lists. We also 
experimen ted with other source ex tractors such as SEX- 
tractor l|Bertin fe Arnoutill996H . and all gave similar 
uncertainties. 

The astrometric solution is provided by ANMATCH, a 
software routinely produces arcsecond precision astro- 
metric solutions for the HATNet/HATSouth observa- 
tions. ANMATCH first uses the engine of astrometry.net 
(jLang et al.l I2010T) for a low order solution (3rd order 
Simple Imaging Polynomial [SIP] tweak) to obtain an 
initial guess, then fits a 3rd order polynomial on a big¬ 
ger matched list between the extracted source and the 
catalogue source to obtain the final solution. 

The histogram of the residuals from the astrometric 
solution, for all 38 CCDs in the K2 Campaign 1 FFIs, 
is plotted in black in Figure [T) The astrometric so¬ 
lution residual is defined as the distance between the 
projected pixel coordinates of catalogue sources and the 
corresponding coordinates for the same stars from our 
source extractor. The median of the astrometric solu¬ 
tion residual for these raw frames is ~ 0.032 pixel for K2 
Campaign 1. 0 Given the Kepler CCD has a plate scale 

5 We expect a factor of two difference between the estimated 
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Fig. 1.— The residual of the astrometric solution for all chan¬ 
nels of the K2 Campaign 1 field FFIs (black) and SFFIs (red and 
hatched). This is computed by comparing the distance between 
the projected catalogue coordinates and the detected source coor¬ 
dinates on the CCDs. 

of 3.98", our astrometric solution residual corresponds to 
0.127". 


2.4. Astrometry for All the Stamps 

The SFFIs are very sparse. For a single channel, typ¬ 
ically more than 95% of the pixels are not downloaded. 
It is impossible to solve for the astrometric solution of 
these SFFIs via a direct catalogue matching. We use the 
FFI astrometric solution as an initial guess to overcome 
this problem. 

We first use FISTAR to extract the sources from the 
SFFIs. We then project the UCAC4 catalogue on to the 
SFFIs with the astrometric solutions obtained in £12.31 
We use an iterative point matching algorithm allowing 
a field centre shift from the FFI to SFFI to match the 
extracted sources and the projected coordinates of cat¬ 
alogue stars. We solve for the distortion between these 
matched pairs using a second order polynomial to obtain 
the final solution. In Figure [2J we show the correspond¬ 
ing region of FFI and SFFI from the same CCD channel 
( K2 Campaign 1, module 13, channel 41). This region 
consists of three stamps in the SFFI observations. We 
marked out the detected source by cyan circles, and the 
projected sources from catalogue by red circles. The orig¬ 
inal K2 targets are marked out in the black circles. Some 
stamps consist of multiple stars. We also show that in 
the top rightmost stamp in Figure [2] the projected cata¬ 
logue indicates that there are additional sources blended 
in the primary source’s PSF, which was originally missed 
by the source extractor, and the light of which would be 
measured together with that of the primary source. 

We compute our astrometric residuals as per £12.31 The 
astrometric residuals on the SFFIs for K2 Campaign 1 
are shown in red in Figure jb The median astrometric 
residual is around 0.034 pixels (0.135"), comparable with 

uncertainties from the source extractor, and our astrometric resid¬ 
ual, given the different methods by which these two uncertainties 
are calculated. We estimated the uncertainty of the source extrac¬ 
tor by compare the .rms. difference of the source positions on two 
frames (allowing a spacial transformation), while the uncertainty 
of the astrometric residual is estimated by the median of residual 
between the extracted position and the projected solution. 


what we achieved on the FFIs. 


2.5. A revised K2 target catalogue 

We projected the UCAC4 catalogue on the K2 Cam¬ 
paign 1 SFFIs using the astrometric solution we obtained 
in 8 12.41 Stars with centroids within 3 pixels from the 
stamp edges were excluded. We also included those stars 
in the original Ecliptic Plane Input Catalogue (EPIC, 
iHuber fc Brvsonl (|2015l P but not in the UCAC4 cata¬ 
logue. The EPIC catal og ue is a com bination of the Hip- 
parcos cata logue (Ivan Leeuwenll2007h . Tyclro-2 catalogue 
ilHog et aLll200Q[h UCA C4 catalogue, 2MASS, and SDSS 
DR9 ( Ahn et aLl I2012T) for the selected K2 target stars. 
There may be systematic offsets between the coordinates 
from the above catalogue, but they are relatively small 
(~ 10 mas), and can be ignored when combining these 
catalogues. We estimated the B and V magnitudes of 
stars not in the EPIC catalogue as per the Kepler In¬ 
strument Handbook. Altogether, we found 14778 stars 
from the UCAC4 catalogue, and an additional 7939 stars 
from the EPIC catalogue only, in K2 Campaign 1. This 
combined set of K2 target catalogue (22717 stars in total) 
is larger by 5% than the total in the original EPIC cat¬ 
alogue (21647 stars in total). This increase will be more 
pronounced for other, more crowded K2 fields. Part of 
the final K2 Campaign 1 target list catalogue is shown 
in Table Q] We provide the centroid positions on the 
corresponding postage stamp, for each target. 


2.6. The refined motion for each module 

VA14 pointed out that the K2 photometry is strongly 
correlated with the centroid positions of the stars. They 
also found that the centroid position of individual stars, 
as determined by their weighted light centres, are often 
not good enough. As such, they chose the centroid mo¬ 
tion of a bright star to represent all the stars observed in 
the same campaign. Taking advantage of our derived as¬ 
trometric solution, we find that by combining many stars 
observed on the same module, we can achieve even better 
constrained X,Y motion tracks. We define the A, Y mo¬ 
tion derived for the centre pixel position of SFFI modules 
as the refined motion for each module. As an example, we 
show in Figure [3] the relative X centroid drift of module 
4, a module in the corner of the focal plane. We did not 
derive a rolling motion for the entire spacecraft to avoid 
correcting for additional rotations between modules. We 
notice that although the spacecraft attempted to correct 
its roll drift every 12 hours, the drifting segments can 
last longer. The drifting segments are defined as a time 
series of smooth X , Y motion without significant out¬ 
liers. We identify each drifting segment, and the outliers 
in between segments, by applying a 1-d edge detection 
method (Sobel operator) on the X motion of each mod¬ 
ule. A example of the edge detection is show for module 
4, in Figure Hi with the red dashed lines separating each 
segment. 


3. PHOTOMETRY 

For each target, we use fiphot d Pa.ll 12012 ) to extract 
photometry in 36 circular apertures around the derived 
centroids. The flux from the sources are estimated by 
summing up all the pixels within an aperture and weight¬ 
ing edge pixels by the fraction of which lie within the 
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Fig. 2. — An 85 X 56 pixel 2 region of the FFI frame (left) and its corresponding SFFI frame (right). The detected sources are marked 
by cyan circles, and the projected sources from UCAC4 catalogue by red circles. The original K2 targets are marked by black big circles 
centred on centroids as determined from our astrometric solution. The white region in the SFFI image were not observed, and are masked 
out. 



Fig. 3.— A segment of relative X centroid drift from module 4, 
between cadence 1000 and 1500. The dashed red lines separate the 
drifting segments we identified via edge detection. 


circular aperture. The background flux is measured by 
taking the median, with iterative outlier rejection, in an 
annulus of pixels around the aperture, then multiplying 
it by the area of the aperture, and subtracting it from the 
flux. The aperture sizes range from 2.5 — 5 pixels, and 
are chosen so as to optimize the photometric precision 
for a wide range of magnitudes. For apertures with sizes 
smaller (larger) than 4.5 pixels, the background annulus 
has inner radii of 5 ( 6 ) pixels and outer radii of 11 (12) 
pixels. They are designed to optimize the photometric 
precision for stars in different magnitude bins. We note 
that for saturated stars (KepMag < 10), our photomet¬ 
ric method cannot capture all the leaked electrons in the 
stamps, therefore leading to degraded photometry. For 
the saturated stars, the fixed aperture approach taken 
by VA14 remains the best way to extract optimized pho¬ 
tometry for now. 

4. LIGHT CURVES AND DETRENDING 

We present the raw aperture photometry light curves 
(described in section Q , and apply a three-step detrend¬ 
ing process on our light curves. The detrending method¬ 
ology is adapted from the HATNet pipeline, as well as 
the Kepler li ght curve detrending pipeline described in 
iHuang et al.1 (|2013l ) . Each step of detrending is aimed to 


correct different aspects of the noise in the light curves. 
Users are able to query light curves detrended up to an 
intermediate step to suit their own purpose 0- In this 
section, we will first describe the properties of our de¬ 
trending methods, and then demonstrate the light curve 
products from each detrending step. At the end of this 
section, we will compare our light curves with those from 
other works. 

4.1. Detrending 

We refer to our light curves from the aperture pho¬ 
tometry as RAWLC. Our detrending pipeline applied to 
these RAWLC can be divided into the following three 
steps: 

(1) External Parameter Decorrelation (EPD); 

(2) Trend filtering (TFA); 

(3) Cosine filtering (COS). 

To correct for the photometric variations due to the 
motion of the spacecra f t, we performed EPD on the 
RAWLC dBakos et al.l 12010) . We follow a similar 
methodology as described in VA14 to deal with the thrust 
fire events of the spacecraft. Instead of correcting for the 
drift effect due to a 6 hour roll, we make use of the drift¬ 
ing pattern of each module we identified from 112.61 We 
first reject the data points that fall in between any drift¬ 
ing segments. We then divide the data into 6 segments 
before detrending, as defined in Tabled These segments 
are designed to separate large amplitude flux offset in the 
data, and allowing each segments to be represented by 
low order smooth functions. For each segment, we itera¬ 
tively fit a 3rd order B-spline through the median mag¬ 
nitude of each drifting segment with 3-er outlier rejection 
until the fit converges. This long term trend represented 
by the B-spline is then removed. 

We then fit for the variation due to spacecraft drift as 
per the following: 

f(m) = Co + Ci sin(27r A") + C 2 cos(27r X) 

+ C 3 sin(27r Y) + C 4 cos(27r Y ) 

+ c 5 sin( 47 r A') + C 6 cos( 47 r A) 

+ C 7 sin( 47 rF) + cs cos(4wY), 

6 http://k2.hatsurveys.org/ 
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in which, X , Y represent the relative X, Y drift of the 
module on which the target sits. The fitted X , Y trend 
is then removed from the original RAWLC, and the 13- 
spline long term trend added back in. This preserves the 
long-term trend while minimizing the effects of short¬ 
term spacecraft motion. The light curves at this stage is 
called EPDLC. 

The shared systematics between the stars are then cor¬ 
rected using an adaptation of the Trend Filtering Algo¬ 
rithm (TFA) designed for Kepler. The idea of TFA is to 
select a set of template light curves, that is representa¬ 
tive of all the systematic variations present in the data. 
Each target light curve is then corrected based on a lin¬ 
ear filter that identifies the shared trends between the 
target and the template light curves. We found that us¬ 
ing only template stars observed in the same channel as 
the target provided the best results. Since the number of 
stars observed in each channel in K2 Campaign 1 is quite 
small, we use all the stars but the target as templates in 
the TFA procedure. The TFA filtered light curves are 
denoted as TFALC. 

The last step is to filter all the low frequency vari¬ 
abilities (mostly due to intrinsic stellar variability) us¬ 
ing a set of cosine and sine funct ions. This method 
was implemented by iHuang et al.1 (120131 1 for the inde¬ 
pendent search of planetary candidates in the original 
Kepler data. We aim to keep all periodicities at or below 
the protected timescale of the transit undistu rbed, while 
minim izing any other variations following Kippin g et al.l 
(|2013T) . The cosine function detrending is applied to the 
EPDLC light curves, the resulting light curves are called 
the COSLC. The cosine function detrending process is 
independent of the TFA process above. Due to its pur¬ 
pose, astrophysical signals such as stellar pulsations are 
no longer preserved in the COSLC. 

4.2. Light Curve Products 

We provide two types of measurements about the pre¬ 
cision of our light curves. We use the point to point 
median scatter around the median (MAD) to represent 
the overall variability in the light curves which is used 
as an estimation of noise level in our transit search al¬ 
gorithm. We also report the 6.5 hour precision as per 
VA14, which characterizes the noise of light curves at a 
time scale relevant to the transit duration of an earth 
analog. Figure [I] shows the MAD of our light curves. 
The best precision of RAWLC, EPDLC, COSLC and 
TFALC for the bright stars are 1.2 x 10 -4 (120 ppm), 
~ 6x 10~ 5 (60 ppm), ~ 5x 10 -5 (50 ppm) and ~ 5x 10~ 5 
(50 ppm), respectively. Figure [5] shows the estimated 6.5 
hour precision of our light curves. The best precision of 
RAWLC for the bright stars are ~ 2 x 10~ 5 (20 ppm), 
and ~ 1.5 x 10 -5 (15 ppm) for all the other three types 
of light curves. We also overlaid the estimation of the 
bottom e nvelo p of the o riginal Kepler 6.5 Hour precision 
based on I Jenkins et al.l (120101 1. 

We show that the EPD process always improves both 
the short time scale (6.5 hours) and the long time scale 
(whole campaign) precision compared to the RAWLC. 
We find a greater improvement for the bright stars com¬ 
pared to the faint stars. The COS filtering process im¬ 
proves the precision in both time scales compared to 
EPDLC. For most stars, the TFALC has a similar or 


marginally worse precision compared with the COSLC, 
but the TFA process tends to preserve the intrinsic vari¬ 
abilities of the stars. We note that for a small fraction of 
the stars, the COSLC and TFALC can have a worse 6.5 
hour precision compared to the EPDLC. This is because 
the cosine filter and TFA algorithm both use linear least 
square method aiming to minimize the overall point-to- 
point scattering in the light curves, which is sometime 
achieved at the cost of increasing noise at specific time 
scale. We compare the noise properties of the light curves 
in Figure [6] by showing the ratio of per point root mean 
square (RMS) and MAD versus magnitude. If the noise 
is composited with pure white noise, this ratio should be 
\J2. The TFALC have the most white noise compare to 
other detrending stages. 

5. COMPARISON WITH OTHER WORKS 

Many teams have developed methods to improve K2 
photometry. We summarize the different approaches ac¬ 
cording to their photometry and detrending methods. 

5.1. Photometry Methods of previous works 

All the teams use aperture photometry method to ex¬ 
tract the light curves from K2 data. However, they differ 
in the details of aperture choice and centroid measure¬ 
ments: 

• Fixed Mask method: 

The fixed mask method is such that the flux of 
the target is summed up over pixels within a 
fixed predetermined mask, while the pixels are ac¬ 
counted in a binary way. VA14 used a combina¬ 
tion of approximate circular aperture and fitted 
apert ures usin g the Keple r Pixel Respo nding Func¬ 
tion (Brvso n et al.| |2010f). | Foreman-Mackev et al.l 
(12015 ) and lAngus et al.l (I2015T) used approximate 
circular but binary apertures (do not include pix¬ 
els partially) and p resent the_photometry from 
the best apertures. iLund et alf (j2015H used the 
density-based spatial clustering of appli cation s 
with noise routine (DBSCAN, lEster et al.1 (|1996l) l 
to chose their pixel mask (aperture). 

• Moving Circular Aperture method: 

lAigrain et al.1 (1201511 used 6 circular apertures to 
extract photometry. The apertures are soft-edged, 
in the sense that pixels straddling the edge of the 
aperture contribute partially to the flux. 

• Centroids from astrometric solution: 

lAigrain et al.1 (|2015h used centroids derived from 
their own astrometric solution with the 2MASS all¬ 
sky point-source catalogue. 

• Centroids from WCS header: 

iForeman-Mackev et al.l (1201511 and lAngus et al.l 
( 201511 used centroids from the WCS header of the 
K2 target pixel files. Only one WCS solution is 
given for the entire time series of each star. 

• Centroids from weighted centre of flux: 

VA14 and ILund et aU (1201511 used the weighted 
centre of flux as the centroids of the stars. We 
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Fig. 4.— The point to point median standard deviation around the median (MAD) versus Kepler magnitude of all the light curve 
products at different detrending stages. From left to right, we show the RAWLC the EPDLC, the COSLC the TFALC and the VALC. 
The dashed red line is the fitted function for the magnitude versus the median MAD in the magnitude bin. The solid horizontal line 
indicates a scatter of 10“ 4 (100 ppm). The vertical scale is logarithmic, and is the same for each panel. 



note, in the subsequent detrending, VA14 used the 
centroid of star EPIC 201611708 instead of the cen¬ 
troids of individual stars. 

In this work we used 36 moving circular apertures, with 
the centroids of apertures determined by precise astro- 
metric solutions, to determine the photometry of each 
star. 

5.2. Detrending Methods of previous work 

There are three different types of “detrending” meth¬ 
ods used by other authors. 


lAigrain et all (120151 ) and iCrossfield et al.l (|2015[ ) 
used Gaussian process model, with the rolling an¬ 
gle as the input variable to detrend the light curves. 
They assume the systematics can be modelled as a 
function form of the rolling angle, and that func¬ 
tion’s form can vary from star to star. 

• Not Detrending: 

iForeman-Mackev et al.l (|2015l) and lAngus et all 
(120151) choose to not detrend their light curve prior 
to the search of signal, but instead, they simul¬ 
taneously fit for the systematics and the signal of 
interest. 


• Decorrelation: 

VA14 , lArmstrong et all (120141) and iLund et all 
(2015) used a self-flat-field method to decorrelate 
the aperture photometry from centroid position 
of the image. There are, however, subtle differ¬ 
ences between these studies. VA1 4 used the cen¬ 
troids from a representative star, lArmstrong et ahl 
(2014J ) seemed to use cen troids for individual stars, 
while ILund et all (120151 ) used the weighted light 
centroids derived for individual stars. VA14 used 
a 1-d decorrelation al o ng th e trajectory of the 
drift, lArmstrong et al.l (|2014D used 2 -d centroid 
surfac e to decorrelate with the flux, and lLund et all 
(120151) used both the 1-d and 2-d approach in their 
pipeline. 

• Gaussian Process: 


In this work, we applied three stages of detrending. In 
the first stage, we applied a similar method as the decor¬ 
relation detrend in VA14. We additionally applied TFA 
and COS filtering to further filter the data. TFA is aimed 
to correct for shared systematics between the stars ob¬ 
served on the same channel, while preserving the stellar 
variability. The COS filtering method aimed to correct 
for any variability in the light curves and is optimized 
for searching for transit signals. 

5.3. Centroids Determination 

We took a similar approach as lAigrain et all (2 0151 ) in 
the determination of centroids by deriving an astrometric 
solution for each image. We made use of a more precise 
catalogue, UCAC4 instead of 2MASS, and the Full Frame 
Image as a better initial guess, and achi eved higher pre¬ 
cision in our astrometric solution. lAigrain et al.1 (|2015l ) 
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Fig. 5.— The 6.5 hour precision versus Kepler magnitude of all the light curve products at different detrending stages. From left to 
right, we show the RAWLC the EPDLC, the COSLC the TFALC and the VALC. The solid black horizontal line indicates a scatter of 
10 —4 (100 ppm). The vertical scale are the same for each panel, linear and in units of parts per million (ppm). The dashed red line is the 
fitted function for the magnitude versus the median 6.5 hour precision in the magn itude bin. The d ashed blue line is the fitted function 
indicate the bottom envelope of the original Kepler 6.5 Hour precision based on I Jenkins et all (120100 . 
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Fig. 6.— Ratio of point to point RMS and MAD versus Kepler magnitude. From left to right, we show the RAWLC the EPDLC, 
the COSLC and the TFALC. The red horizontal line indicates the value of v2, which should be the value or their ratio for pure white 
noise. 


reported a typical root mean square of the astrometric 
solution of 0.4", or approximately 0.1 pixel, ~ 3 times 
larger than our typical astrometric residuals. 

5.4. Photometric Precision 

To date, only VA14 have released their detrended K2 
Campaign 1 light curves, therefore we will focus on com¬ 
paring our photometry precision with their work. 

We show in Figure [7] the precision ratio between our 
light curves and the VA14 light curves, for the same 


stars at both times scales. The EPDLC, COSLC and 
TFALC from this work have comparable precision com¬ 
pare to VA14 light curves on the 6.5 hour time scale, and 
a smaller point-to-point scatter over the entire observa¬ 
tion length. 

lAigrain et al.l (|2015f l presented their ctmad (similar to 
MAD) and 6.5 hour Combined Differential Photometric 
Precision (CDPP) for K2 engineering data photometry. 
Their best precision for ctmad is ~ 300 ppm for the bright 
stars, and 60 ppm for 6.5 hour CDPP. In this work, we 



























Kep Mag Kep Mag Kep Mag Kep Mag 

Fig. 7.— Photometry precision of our light curves compared to VA14 K2 Campaign 1 light curves. Top panel: 6.5 hour precision ratio 
between our light curves and VA14 light curves versus Kepler magnitude. Bottom panel: Per point MAD ratio between our light curves 
and VA14 light curves versus Kepler magnitude. From left to right, we show the RAWLC the EPDLC, the COSLC and the TFALC. 
The red horizontal line indicates the value of 1. We note the vertical scale in the two panels are different, and both in log scale. 


achieved higher precision on both time scales (50 ppm 
and 15 ppm, respectively). Although, we caution that 
the noise characteristic for K2 engineering data and K2 
Campaign 1 data could be different. 

5.5. Power Spectrum 

Previous works dLund et al.l 1201.51 lAngus et, all 1201511 
noted the detrended light curves from a ID decorrela¬ 
tion may still have residual spikes around the harmonics 
of ~ 47.2271 pKz in their power spectra. These residu¬ 
als may be largely due to aliasing of the low frequency 
power, induced by data gaps from rejected points during 
thruster firing. These harmonics are less obvious in our 
light curves. 

We computed the Discrete Fourier Transfor- 
mation (DFT) pow er spectrum using vartools 
(jHartman et_aLj _[2008jl for t he star EPIC 201183188, 
following I Angus et al.l (|2015H . We find our EPDLC 
does not show the frequencies corresponding to the 6 
hour roll in the DFT power spectrum present in the 
VA14 light curves. Figure [5] compares our EPDLC for 
EPIC 201183188 with the VA14 reduction, and Figure 
[9] compares the DFT power spectrum of our EPDLC 
with that of the VA14 light curve. The offset around 
epoch BJD 2456015 in the VALC could be blamed for 
contributing to noise peaks corresponding to the space 
craft rolling frequency. 

We further investigate the noise characteristic of the 
light curves by computing the median of DFT power 


spectra of 1661 stars in the magnitude range of 10-12. 
To eliminate the influence from the long term trend, we 
first filter out the strongest low frequency peak in the 
light curves, and then recompute the DFT power spec¬ 
tra of each star before taking the median. We show this 
median spectrum in Figure [TO] The peaks related to the 
rolling frequency are still present, but have less power 
compared to the median DFT spectrum computed with 
the same set of stars from VA14. When the signal from 
stellar variability is strong, such as in the case of EPIC 
201183188, the systematic noise peaks are negligible. 

5.6. Example Light Curves 

To demonstrate the products of our photometric ex¬ 
traction and detrending in the context of stellar variabil¬ 
ity and transit searches, we show example light curves 
for stars with known variability and transiting signals. 

EPIC 201711881, is a Cepheid c andidate discovered 
by the AS AS project dSchmidt et al.l[2009h . The original 
discovery paper reported a period of 2.7353±5 x 10~ 4 
days. The period we detected is ~ 2.735 day, consistent 
with dSchmidt, et aDl2009ll . The RAWLC and EPDLC 
folded with the Cepheid period, are shown in Figure fill 
Our RAWLC without any detrending, already shows a 
clean periodic signal. In addition, our EPDLC preserves 
the amplitude of the pulsation after removing the sys- 
tematics. We also find a periodic eclipsing signal in the 
light curve with twice the pulsation period. This eclips¬ 
ing signal can be visually seen, even in the RAWLC. 
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Fig. 8.— The EPDLC for EPIC 201183188 from this work 
(black), and VA14 (red). The DFT power spectrum of this star is 
shown in Figure [9] 
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Fig. 9.— The DFT power spectrum for EPIC 201183188 from 
this work (top panel) and VA14 (bottom panel). Note the vertical 
scales are different in the two panels. The dash lines indicate the 
corresponding frequency associated wi th the spacecraft roll motion. 
For more discussion, see I Angus et a ni f2om 



Fig. 10. The median DFT power spectrum of 1661 stars in 
the magnitude range 10-12. Red dashed lines show the result from 
the VA14 light curves. Black solid lines show the result from the 
EPDLC of this work. In the small window of this figure, we show 
the zoom-in of these spectra around the space craft roll frequency 
~ 47.2271p, Hz. 


The COSLC, phase folded with the eclipse period, is 
also plotted in the bottom panel of Figure [IT] It seems 
that this system is more likely to be an eclipsing binary 
system with large out-of-transit variation rather than a 
real Cepheid. To further characterize this signal, the 
COSLC from the general detrending pipeline is not good 
enough. The COSLC we show in Figure fill have been 
reconstructed after the discovery of the eclipsing signal. 
Only the out-of-eclipse part have been filtered by the co¬ 
sine filters in order to preserve the shape and amplitude 
of the eclipsing signal. 

WASP-85b: WASP-85 (EPIC 201852715), was ob¬ 
served in K2 Campaign 1, in module 15, channel 49. 
1182 other stars were observed in the same module. We 
show in Figure [12] the light curves of WASP-85 at dif¬ 
ferent detrending stages from this work. In the second 
panel from the top, we also overlaid the VA14 detrended 
light curve. WASP-85b has a period of ~2.65 days, and 
known depth of ~ 1.6%. We show the WASP-85 light 
curve folded with the detected period and epoch in phase 
space for both the COSLC and TFALC in the bottom 
of Figure [TU 

K2-3: K2-3 (EPIC 201367065), was observed in K2 
Campaign 1, in module 12, channel 40. The host star 
is an M dwarf, with three tra nsiting super-earths discov¬ 
ered bv lCrossheld et all (1201511 . We show the light curves 
for K2-3 in Figure H3l The transits of the biggest planet 
(1 mrnag) is visible in our RAWLC, and the transits of 
all three planets are visible in all the other light curves. 
We also show the phase folded COSLC and TFALC for 
all three planets in the bottom panel of Figure [T51 

EPIC 201613023: This star was identified as a tran¬ 
siting planet candidate system bv lForeman-Mackev et al.1 
(120151) . We show our light curves in Figure [Si The 
transit signal has depth of 400 ppm. Individual transits 
from the planets are visible in the EPDLC. TFALC and 
COSLC. The phase folded COSLC and TFALC with 
the detected epoch and period are shown in the bottom 
panel. 


6. CONCLUSION AND DISCUSSION 

In this article, we present our effort to extract high 
precision photometry from K2 Campaign 1 data. Our 
method has three distinct advantages: 

• Making use of accurate astrometric solution 
(0.127" or 0.034 pixels) from the FFIs for aperture 
centroiding; 

• Providing photometry for all sources on the 
stamps, not only for the proposed targets from the 
input catalogue; 

• Presenting light curves with very low systematic 
variations. 

Our extracted light curves are of high precision at both 
the long (entire campaign) and short (6.5 hours) time 
scales, even for the raw light curves without any detrend¬ 
ing. Light curves derived from all 36 photometric aper¬ 
tures at all four detrending stages are provided for the 
public at http : //k2 . hatsurveys . org 
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Fig. 13.— Same as Figure [T2l but for K2-3 (ICrossfield et alii2019) . In the bottom panel, we show all the three planets detected in this 
system folded with their own periods and epochs. K2-3 b (P~ 10 day, R~2 R® ), c (P<~ 25 day, H~1.7Rq ), d (P<~ 44 day, R~1.5 R® ) 
are presented in red,blue and black curves, respectively. We note that the vertical scales in different panels are different. 



































A Mag A Mag A Mag A Mag A Mag 


13 


EPIC 201613023 


12.045 [- [Raw LC 
12.047 
12.049 
12.051 


12.045 

12.047 

12.049 

12.051 


12.048 - 

12.05 
12.048 

12.049 

12.05 



EPD LC 
VA LC 


ITFA LC 


» /# M | 


(A) 




(B) 




(C) 







1970 1980 1990 2000 2010 2020 2030 2040 2050 2060 

BJD-2454833 (Day) 



bo 

cd 



Hours from Mid-Transit 


Hours from Mid-Transit 


Fig. 14.— Same as Figure fl2] but for EPIC 201613023 (|Foreman-Mackev et al.|[2015D . In the bottom panel we show the detected planet 
candidate folded with its period (~8.28 days). The depth of the transit is ~ 400 ppm. We note that the vertical scales in different panels 
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TABLE 1 

K2 TARGET LIST 


UCAC4ID 

RA 

DEC 

J 

H 

K 

B 

V 

g 

r 

i 

x(to) b 

y(to) b 

channel c 

K2ID d 

flag 3 

UCAC4-555-033290 

101.764743 

20.947725 

11.705 

11.192 

11.067 

11.927 

11.661 

11.756 

11.696 

11.726 

41.098906 

885.251787 

24 

202071861 

AAeeeAAAe 

UCAC4-555-033327 

101.789413 

20.952493 

14.487 

12.700 

12.157 

15.417 

14.481 

14.920 

14.144 

13.857 

42.030470 

863.957681 

24 

202071849 

AAeeeAAAe 

UCAC4-555-033328 

101.789718 

20.941804 

15.207 

14.082 

13.808 

15.967 

15.332 

15.603 

15.144 

14.984 

32.436217 

865.242115 

24 

202071849 

eeeeeeeee 

UCAC4-555-033330 

101.790216 

20.942716 

16.218 

14.420 

13.828 

17.923 

16.773 

17.021 

16.130 

15.738 

33.183423 

864.694662 

24 

202071849 

eeeeeeeee 

UCAC4-555-033335 

101.795353 

20.945481 

11.554 

11.266 

11.224 

11.589 

11.466 

11.478 

11.559 

11.677 

34.961871 

860.006732 

24 

202071849 

AAeeeAAAe 

UCAC4-555-033336 

101.795597 

20.950077 

16.008 

14.560 

14.051 

17.444 

16.468 

16.771 

16.030 

15.720 

39.036805 

859.140643 

24 

202071849 

eeeeeeeee 

UCAC4-555-033338 

101.799223 

20.951700 

15.611 

14.331 

14.073 

15.935 

15.416 

15.625 

15.249 

15.125 

39.999106 

855.879192 

24 

202071849 

eeeeeeeee 

UCAC4-555-033418 

101.876901 

20.980552 

14.091 

13.383 

13.252 

14.427 

14.080 

14.193 

14.071 

14.044 

55.344496 

786.888096 

24 

202068459 

AAeeeAAAe 

UCAC4-555-033420 

101.876908 

20.972664 

12.859 

11.979 

11.702 

13.485 

12.948 

13.181 

12.839 

12.712 

48.289514 

788.011635 

24 

202068459 

AAeeeAAAe 

UCAC4-555-033425 

101.881285 

20.975126 

16.224 

13.317 

12.709 

17.317 

16.082 

16.193 

15.201 

14.762 

49.903488 

784.005099 

24 

202068459 

eeeeeeeee 


a The photometry flag indicator: A-APASS photometry; T-Tycho photometry; e-estimated with 2MASS photometry 

b the X and Y coordinate of star at time 0, time 0 is defined as cadence 1 (BJD-2455895.528) for K2.0, and cadence 102 (BJD-2455975.178) for K2.1 
c the channel number of which the star is observed 
d the given K2 ID of the star 
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TABLE 2 

Light Curve Segments used in EPD 


Segment No 

start Cadence 

end Cadence 

i 

0 

454 

2 

455 

1005 

3 

1006 

1989 

4 

2050 

2314 

5 

2315 

2997 

6 

2998 

4020 



